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1 Introduction 

The analysis of the structure and properties, as well as the evaluation of exact form 
factors, is one of the central problems in 1+1 dimensional quantum field theories. 
One of the main reasons for their distinct role is that they serve to compute very 
efficiently correlations functions of local operators 0{x). Instead of a perturbative 
expansion in the coupling constant one may expand the correlation functions in 
terms of exact expressions of n-particle form factors, that is the matrix element of 
a local operator 0{x) located at the origin between a multiparticle in-state and the 
vacuum 

The operators Zj^{9) are creation operators for a particle of type /x as a function of 
the rapidity 9. 

Since the original proposal of this method to evaluate correlation functions fl]], 
various schemes have been suggested to compute these objects. One of the origi- 
nal approaches is modeled in spirit closely on the set up for the determination of 
exact scattering matrices. It consists of solving a system of consistency equations 
which have to hold for the n-particle form factors based on some natural physical 
assumptions, like unitarity, crossing and bootstrap fusing properties |l|, |^ 

^p'^\■■■p^i^^i+l■■■ f a a \ — Tp^\---t^i+ii^i--- / o o \q (q \ (o\ 

^o|M,.../.„(^^ + 27rz,...,^„) = 7?,i^r''^^-^"'^n^2,...,^n,^i), (3) 

n 

_Res F;}^r'-^-{e + m, e,, 6,. . . 9,,) = ^(1 - 7? J] S,,,{eoi))F^^^-^-{e,. . . 9n). (4) 

Here 7^ is the factor of so-called local commutativity defined through the equal 
time exchange relation of the local operator 0{x) and the field O^ijj) associated to 
the particle creation operators Zjj.{0) 

0,{x)0{y)=^''^0{y)0,{x) for x^ > y\ (5) 

The factor 7^ is very often omitted in the analysis or simply taken to be one, but it 
can be seen that already in the Ising model it is needed to set up the equations con- 
sistently 0] . A consequence of its presence is that a frequently made statement has 
to be revised, namely, that (H)-® constitute operator independent equations, which 
require as the only input the S-matrices SijiOij) between particles of type i and j 
as a function of the rapidity difference 6ij = 6*, — Oj. In the following manuscript we 
demonstrate that apart from ±1, which already occur in the literature, this factor 
can be a non-trivial phase. Thus the form factor consistency equations contain also 
explicitly non-trivial properties of the operators. 



Trying to find solutions to these equations has been pursuit successfully for many 
models and has led to the determination of closed exact expressions for n-particle 
form factors for a wide class of local operators 0{x), e.g. |]l], Q. 

Alternatively some authors develop methods which borrow ideas which have 
proven to be very powerful in the context of conformal field theory, where the use 
of symmetries and their related algebras has led to a successful determination of 
correlation functions 0]. Yet, the most direct way to compute the matrix elements 
in (|l]) is to find explicit representations for the operators ZU6) and 0{x). For 
instance in the context of lattice models this is a rather familiar situation and one 
knows how to compute matrix elements of the type ([l|) directly. The problem is then 
reduced to a purely computational task (albeit non-trivial), which may, for instance, 
be solved by well-known techniques of algebraic Bethe ansatz type, e.g. 0. In the 
context of field theory a similar way of attack to the problem has been followed by 
exploiting a free field representation for the operators ZU6) and 0{x), in form of 
Heisenberg algebras or their q-deformed version. So far a successful computation of 
the n-particle form factors with this approach is limited to a rather restricted set of 
models and in particular for the sine-Gordon model, which is a model extensively 
studied by means of other approaches [^ ^, only the free Fermion point can be 
treated successfully [^ |[ so far. One of the main purpose of this manuscript is to 
advocate yet another approach, namely the evaluation of the matrix elements (|l|) 
based on an expansion of the operators in the conventional fermionic Fock space. 
Recalling the well-known fact that in 1+1 space-time dimensions the notions of 
spin and statistics are not intrinsic, it is clear that both approaches are legitimate. 
Since the model we mainly consider in this manuscript, the Federbush model, is 
closely related to complex free Fermions the usage of fermionic Fock operators 
seems natural. Nonetheless, we expect this procedure to hold in more generality 
and to allow an extension to other models. 

Our manuscript is organized as follows: In section 2 we recall an explicit fermionic 
free field representation for the particle creation operators Zj^{6) occurring in (|l|) 
valid for all diagonal scattering matrices. In section 3 we treat the complex free 
Fermion, we provide a generic expression for a potentially local operator and spe- 
cialize it to particular operators whose form factors we directly compute, namely the 
order and disorder field and various components of the energy-momentum tensor. 
In section 4 we extend this analysis to the Federbush model and show in partic- 
ular how it is related to homogeneous sine-Gordon (HSG) models on the level of 
the scattering matrix. In addition we analyze the momentum space cluster prop- 
erty. We pay special attention to the factor of local commutativity. In section 5 
we propose a Lie algebraic generalization of the Federbush model. In section 6 we 
sustain the relation between the Federbush and the HSG-models by carrying out 
the ultraviolet limit. We state our conclusions in section 7. 



2 Fock space representation for the FZ-operators 



In order to proceed in the way as outlined above, we have to provide exphcit rep- 
resentations for the creation operators ZU6) and the fields 0{x). The former oper- 
ators are characterized by their braiding behaviour, i.e. when they are exchanged 
they pick up the scattering matrix as a structure constant. We restrict our consid- 
erations in this manuscript to theories in which backscattering is absent, such that 
the exchange algebra for the Z-operators reads [§ 



zlmzUe, 



s,,{e,,)z]{e,)z\ 



exp[27ri5,,(%)]ZJ(^,)4( 



(6) 



As indicated in @, the scattering matrix Sij{9ij) can be expressed as a phase. 
Identical relations hold for the annihilation operators, i.e. Z'^{6) -^ Z{6) in (^. 
When we braid a creation and an annihilation operator the presence of an additional 
central term was suggested in [|10| 



h)z]{e,) 



s,,ie,^)z]ie,)z. 



2n6,,6{9, - 9,) , 



(7) 



which ensures that one recovers the usual (fermionic) bosonic (anti)-commutation 
relations in the case (5* = —1) S = 1. The relations (§|)-(0) are commonly referred 
to as Faddeev-Zamolodchikov (FZ) algebra. A representation for these operators in 



the bosonic Fock space was first provided in |TT 



zm 



exp 



I / d9' Suie - 9')alie')aiie') 



4(0) 



(8) 



By replacing a constant phase with the rapidity dependent phase 5ij{6) and turning 
the expression into a convolution with an additional sum over /, the expression (^) 
constitutes a generalization of formulae found in the late seventies |]12[, which in- 



terpolate between bosonic and fermionic Fock spaces for arbitrary spin. The latter 
construction may be viewed as a continuous version of a Jordan- Wigner transfor- 
mation |jl3|, albeit on the lattice the commutation relations are not purely bosonic 



or fermionic, since certain operators anti-commute at the same site but commute 
on different sites. Alternatively, one may also replace the bosonic a's in (H) by 
operators satisfying the usual fermionic anti-commutation relations 



{a0),a,{e')} = O 



and 



a^{e),aUe')} = 27rS^,s{e-e') 



(9) 



and note that the relations (y) are still satisfied |1J]. In the following we want to 
work with this fermionic representation of the FZ-algebra (^). Having obtained 
a fairly simple realization for the Z-operators, we may now seek to represent the 
operator content of the theory in the same space. In general, this is not known and 
we have to resort to a study of explicit models at this stage. 



3 Complex free Fermions 

To illustrate the procedure, to fix some of our notations and to set the scene for 
the Federbush model, we will commence with the free Fermion. Let us consider N 
complex (Dirac) free Fermions described as usual by the Lagrangian density^ 

N 

Cpp = J2 ^ai^Yd^, - m„)Va • (10) 

The associated equations of motion, i.e. the Dirac equations {i'j^d^ — ma)ipa = 0; 
may then of course be solved with the help of the well-known Fourier decomposition 
of the complex free Fermi fields 

^^{x) = I -^ {a.{p)u^{p)e-'P--' + 4(p)t;„(p)e^^-^) a = l,...,N. 

fii) 



We abbreviated as usual \/'ni\ + v\ = Pa ^^^ employed the Weyl spinors 



Mp) = J^ ( ^^e/2 j and v^ip) = iJ^ I ^^e/2 j • (12) 

The amplitudes of the scattering matrices are simply Saa' = —1 for all combinations 
of a, a' and the FZ-algebra coincides by construction with the Clifford algebra (^, 
that is Za{0) = aa{6). 

A further property, which we want to exploit here and in the next section, is 
the [/(l)-symmetry of the Lagrangian £ff, that is changing 

i^ai.^) -^ Va^ai^), (13) 

with rj^ G f/(l) leaves the Lagrangian in (|Ty) invariant. This simple symmetry will 
allow an a priori judgement about vanishing form factors. 



"We use the following conventions throughout the paper: 
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3.1 Form factors of some local operators 

Let us now define a prototype auxiliary field 



1 



K'^ie^e') faj,(0)4(^')e'^^"''''^" + a^{e)aa{e')e-'^^+p'^-'^ 



X^A^) = -^ / dOdO' 

+k"(^, e' - m) (ai{e)a^{e')e'^^~'p'^-'' - a^(^)a1,(^')e"'^^"^'^"^)] • (14) 

This field is essentially bilinear in the free Fermi fields up to the function n{6,6'), 
whose precise expression, which gives the field its individual characteristic, we will 
leave generic for the time being. The properties of this function, like k{6 — in, 6' — 
in) = K,{6,d'), as well as the form of the space-time dependent exponentials, are 
dictated by the crossing in ([l|). It means bringing consistently some of the particle 
creation operators ^^(^) to the left of the operator 0(0) introduces these con- 



straints. Fields of this nature appear already in |T6|. We now want to compute the 
matrix element of a general operator composed out of Xk(^) 

0xS(^) = :e^SW. . (15) 

The direct computation of matrix elements related to these fields is straightforward 
by employing Wick's first theorem^. Noting that the contribution from the normal 
ordered part is of course zero, since all annihilation and creation operators are 
brought to the left and right, respectively, we obtain for instance 



(16) 



F2^"i^"(^i,^2) = -^ I dede'n''ie,e')la^ie)a^ie')aUOi)aii02 
= I dede'K^ie, e')5{o - 92)5(0' - 9i) = ^"(^2, Oi 

Proceeding in this way to higher particle numbers, we compute 

1 r " 

F2r'"^""(^i, . . . , O^n) = - d9[... d9',^ n '^"(^2.-1, ^2.) detV'- , (17) 

where V^ is a rank i matrix whose entries are given by 

V'^^ = cos'[{z-j)7r/2]6{9'^ - 9j) , l<t,j<i. (18) 

We used the identity 



aa{9'J ■ ■ ■ as.{9',)aM)aJ.{0[)al{9i)ai{92)aU03) ■ ■ ■ ai(^„) + 



aa{0'n) ■ ■ ■ as.{9',)a^{9'2)a^{9[)aU9,)ai{92)aU9,) ■ ■ ■ ai{9^) + ■ ■ ■ = Det P". (19) 



^The difference between the product of some linear operators and its normal ordered prod- 
uct has to be a c-number determined by all possible contractions, i.e. for the linear operators 
A,B,C',... holds ABC ... — -.ABC . . . : = sum over all possible contractions (see e.g. [|l9|). 



A further generic field, which we want to study and which, in contrast to O^" , now 
possesses non-vanishing matrix elements with an odd particle number is 

O'^^ix) = :^„(x)e^"(^'): . (20) 

This field involves the fermionic field with the spinor structure stripped off 



i^ai^) 



dpi 

27rpO 



a„(p)e-'P-" + 4(p)e^P-" 



(21) 



Similarly as before we compute the matrix elements 



n\ 



^2r'"^"^'"^(^i,---,^2n+i) = ^ I de[...de',^,^,l[K-{e',,,e',^^,)detv'-^K (22) 



j=i 



Note that 0^^{x) and 0^"{x) are in general non-local operators, in the sense 
that it is not guaranteed that they (anti)-commute for space-like separations, i.e. 
[0{x),0'{y)] = for (x — yY < 0. At the same time F^ is just the matrix ele- 
ment as defined on the r.h.s. of ([1|) and not yet a form factor of a local field, in 
the sense that it satisfies the consistency equations (H)-®, which imply locality of 
00. In order to distinguish between this two different situations we denote matrix 
elements in general by F^ and form factors of local operators by F^. For instance, 
as a consequence of the monodromy equation (^, a necessary condition for these 
two functions to coincide for Xk(^) i^ 



.e' + 2m) = -ji^K''{d,d') 



(23) 



Before specifying the functions k more concretely such that the corresponding O's 
become local, we would like to compare briefly the generic operators of the type 
(|T^), (|TB|) and ( pOD with some general expressions for "local" operators which ap- 



pear in the literature [|T^, |TJ, [T^ . We carry out this argument in generality without 
restriction to a concrete model. Let us restore in equation (|l|) the space-time depen- 
dence, multiply the equation from the left with the bra- vector (Zj^^(9n) ■ ■ ■ Zj^^{9i) 

and introduce the necessary amount of sums and integrals over the complete states 
such that one can identify the identity operator I 



n=l...oo 

E 

n=l...oo 

{0{x) I. 



dOi... de. 



n\{2Ti) 



dOi... dOr. 
n\{2-nY 



-F^\^^-^^^ie,...e^)(zl{e^) 



-zli 



o{x)zl{e,)...zl{e^)){zl{e^)...zl{e,) 



'i'EjPj-^ 



^A rigorous proof of this statement to hold in generaUty is stiU an open issue. 



Cancelling the vacuum in the first and last line, and noting that we can replace 
the product of operators, which is left over also by its normal ordered version, we 
obtain the expression defined originally in ||l7 



n=l...oo ^ ' 



Ml--/^n "°° 



Hence this field is simply an inversion of (|I|). From its very construction it is clear 
that 0{x) is a meaningful field in the weak sense, that is acting on an in-state 
we will recover by construction the form factor related to 0{x). In addition, one 
may also construct the well-known expression of the two-point correlation function 
expanded in terms of form factors, as stated in fl^. However, it is also clear that 



0{x) 7^ 0{x), simply by comparing (|2^ ) and the explicit expressions for some local 
fields occurring in the free fermionic theory, e.g. (0), (^) and (0). The reason is 
that acting on an in-state with the latter expressions the form factors are generated 
in a non-trivial Wick contraction procedure, whereas when doing the same with (|2^) 
the Wick contractions will be trivial. Therefore general statements and conclusions 
drawn from an analysis made on 0{x) should be taken with care. It is also needless 
to say that from a practical point of view the expression (|24|) is rather empty, since 
the expressions of the form factors Fn '^^'"'^"(^0^ . . ,0^^ themselves are usually not 
known and their determination is in general a quite non-trivial task. In |l^, |14|, Q 
the integration in the formula ( P^ is a rather artificial contour integration which 
takes care about analytic continuations of values of in. This does not seem to be a 
fundamental feature, since it remains completely obscure how to incorporate bound 
states in this manner. 

Let us now return to our concrete analysis by specifying k. 

3.1.1 The order and disorder field 

Having in mind to proceed to the Federbush model, we will restrict ourselves from 
now on to the case of two complex Fermions, i.e. iV = 2 in the Lagrangian (|TDp. 
The free fermionic theory possesses some very distinct fields, namely the disorder 
and order fields 

/x„(x) = :e'^"('=): and a^ix) = :V^„(x)/i„(a;):, a = 1,2, (25) 

respectively. The names for these fields result from the ultraviolet limit, see also 
section 5, since then they flow to their equivalent counterparts in the conformal field 
theory [l^, namely to primary fields with scaling dimension 1/16. We introduced 



here the fields 

— — (0—0'^ 

u;^{x)=x:ix), with ^\0,0') = -n^(-0,-0')='--^-^^—-^^ . (26) 



Admittedly, the precise form of the fields u!a{x) appears to be slightly unmotivated 
at this stage. However, we will provide a better rational for this in the next section, 
where we see that they originate by relating a so-called triple normal ordering pro- 
cedure for a field, which can be constructed directly from the Fourier decomposition 
of the free Fermi fields (pA]), to another one associated with the usual Wick normal 
ordering. It will turn out that the field uJa{x) emerges as the limit of a Federbush 
model field to the free fermionic theory, i.e. limA_^i/2 ^^(a;) = ^a{x), see equation 

m- 

Let us now compute the form factors related to the above mentioned fields /iQ,(x) 
and (Jaix). Using the particular form of k°'{9,9') as defined in (p^), we compute 
the integrals in (|l^) and obtain a closed expression for the n-particle form factors 
of the disorder operators 



I'^r-'aniXi, X3, . . . , X2n-l)Bn,n , (27) 



F,T"'\-0u...,-e2n) = (-ifF^';^!-"''!^!,...,^ 



with 



11 (^2i-l ^2i-l) 11 (^ 



2 „2 



_ i<*<i<ri i<i<j<m _ detV"'(x^)detV"(x^) , , 

l<i<j <n+m 

Associated with the particles and anti-particles we introduced here the quantities 
Xi = exp(^j) and Xj = exp(6'j), respectively. The variable Ui can be either of them. 
We also employed the elementary symmetric polynomials ak{xi, . . . , Xn), defined as 

ak{xi,...,Xn) = ^ xi^...xi^ (29) 

Zi<...</fc 

(see e.g. [^ for more properties), the Vandermonde determinant of the rank i 
matrix V^ whose entries are given by 

V.jix) = ixjy-\ l<t,j<i (30) 

and the determinant of the rank i — 1 matrix W^^^ with entries 

W|,-i(x) = a2._,(xi, . . . , X,) , l<t,j<i -I. (31) 



The relations between F^ and F^ as stated in (27) follow most trans- 



parently from ( PBD and (|T^. One may easily verify that the expression ( P7| ) indeed 
satisfies the consistency equations (|D-(^ with 7^" = —1 for a = 1,2. We justify 
this choice in the next section by carrying out the Federbush model —>■ two complex 

8 



free Fermion limit. Noting that Hai^) is invariant with respect to the symmetry 
property ([l3|), it follows immediately that 



^,,\aa...aapp...pp^^^^ . . . ^fc, ^fc+1 • • • , O^^i) =0, « ^ ^, ^ /, «' = 1, 2, 1, 2 . (32) 



This means that, up to a re-ordering of the particles, the expressions reported in ( ]27|) 
are in fact the only non- vanishing form factors related to Hai^) for a G {1,2,1,2}. 
In a similar way we compute the n-particle form factors of the order operator 

pcrjll(nxll)/ n n \ — ( i \n po-2|2{nx22) //] n \ 

^2n+l [-til,. . .,-U2n+l) — [-'-)^2n+l [^^l, ■ ■ ■ , >^2n+l) 

= r2""V„(xi,. . . ,S2n-l)'B„,„+i, (33) 

As a consistency check, one may once again verify that (|33| ) fulfills the form factor 
equations (0)-(|D with 7I" = 1 for a = 1, 2. Again, we postpone the justification of 
this choice to the next section by carrying out the Federbush model — > two complex 
free Fermion limit. Noting that cra{x) -^ Va'^a{x) by (0), it follows immediately 
that 

^a^\aa...a.pp...pp^^^^ • . • , ^., ^.+1 . . . , 9k+l) =0 for « ^ ^, ^ ^ + 1 • (34) 

Of course, this way of proceeding also works for the real free Fermion and one 
may recover the well-known expressions of the literature |TB|, ^, |^ . As a difference 



to our previous computations, however, we have to take care of more contributions 
in the contraction procedure. Keeping the form of k°'{6,6') as defined in (|26|) , but 
taking a = a, we compute for instance 



dOde'K^e, e')[5{e - e2)S{e' - e^) - s{e - ei)S{e' - 62)] 

= i tanh — . (35) 

Proceeding in this way to the higher n-particle form factors we only have to replace 
in (^) the matrix V^ with V^, whose entries are Vfj = 6{6^ ~ Oj). Computing the 
integrals we get 

FL{di, . . . , ^2n) = i" Pf (^) = r\/d^tl4 = «" n *^^^ ~ ' (2^) 

l<i,j<2n 

where A is an anti-symmetric {2n x 2ra)-matrix whose entries are given by Aij = 



tan.h.6ij/2 and Pf denotes its Pfaffian^. In a similar way we compute the n-particle 
form factors of the order operator 

F2Vi(^i,...,^2n+i)=^"PM) = ^" n tanh^. (37) 

l<i,j<2n+l 

Expressions of the type (|36|) and (^) can be found aheady in the first paper of 
|1T6[| . The product expressions for F^ and -F;^+i were also derived in ||2^ and 0, 
respectively, by means of solving the form factor consistency equations (H)-® • 

3.1.2 The energy momentum-tensor 

A further field which plays an important role in any theory is the energy-momentum 
tensor, which for the free Fermion in our normalization simply reads 

n = 2ti:i;,rdu^i : + ■.^2rduM- (38) 



With the help of equation (|ri|) we compute easily 

T^" = xtr + xtr > (39) 

using 
tl^{e, ~e) = -27rima{paT sinh ^^, tl^{9, d) = -2nima{po,Y cosh ^^ , (40) 



where we recall the definition of Xk(^) from (0). A specially distinct role is played 
by the trace of the energy-momentum tensor, since on one hand it is directly pro- 
portional to the operator which breaks the conformal invariance p2| and on the 
other hand it occurs explicitly in various computations associated to the ultraviolet 
limit like the c-theorem p3| and the A-sum rule |2^ (see section 5). It acquires the 



explicit form 

T*"^ = 2imi:^^i/j^: + 2im2:^p2i^2-=Xti + Xt2 ^ t„(6', 6*) = 27rim^ sinh — — . (41) 

It is clear that only the two-particle form factor can be different from zero and we 
compute it in an analogous way as in the previous section, that is using Fourier 



8 Denoting the permutation group of 2n indices by S271 and the signature of the permutation tt 
by sgn^ir), the Pfaffian of a matrix A is defined as 

1 " 



10 



decomposition ([TT]) with subsequent contractions, 

Ff '""(^,^) = -27rim,p'^sinh^^, (42) 

F^^'^^l^,^) = -27rim„p'^cosh^^, (43) 

Fr-'""(^,^) = Ff-^''''{e,d) = -2mmlsmh^^. (44) 

When taking a = a, these expressions coincide with the ones which may be found in 
the hterature for the real Fermion. As usual, we may verify that various equations 
which hold for the operators themselves also hold for the associated form factors. 
For instance, the conservation of the energy-momentum tensor 



d^T'^" 



P^^T^'" 



(45) 



is reflected by the fact that 

[P +P ) F2 = -{P +P )F^ ■ (46) 

Here we used the explicit form of the momentum operator 

P, = r dx'T\ = Y.I ^ iPa), iaiip)a^ip) - a^ip)aUp)) (47) 

J-oo ^^^ J '^'^Pa 

when changing in (|45|) derivatives to commutators by means of the Heisenberg 
equation of motion. It is then easy to verify that [P^,aJ^(p)] = {pa)^i<i\{p) and 
[P^,aa(p)] = —{j>a)iJiO'a{p)i such that we verify exphcitly 

d,x':.{x)=i[P,,xl{x)], (48) 

which is of course what we expect. Equation (^§1) is a further support for the 
consistency of the generic definition of xZ (^) i^ (0) ■ 

4 The Federbush Model 



The Federbush model |25| was proposed forty years ago as a prototype for an 



exactly solvable quantum field theory which obeys the Wightman axioms |^, ^, 



|28[| . Formally it is closely related to the massive Thirring model ||2^. It contains 
two different massive particles \l/i and ^2- A special feature of this model is that the 
related vector currents J^ = ^al^'^a, ct G {1, 2}, whose analogues occur squared in 



11 



the massive Thirring model, enter the Lagrangian density of the Federbush model 
in a parity breaking manner 

Cf=Y, ^aitl^d^ - m,)^, - 2n\s^,J^J^ (49) 

a=l,2 

due to the presence of the Levi-Civita pseudotensor e. It is then easy to verify that 
the related equations of motion 

(lYd^ - rnO^i = 27rA£^, J^V^i, i^^d^ - ^2)^2 = 27TXe,^r,r^2, (50) 
can be solved by 

^1 = ;exp(2v^a02);^i = $^^1, ^2 = iexp(-2v/^a0i);^2 = '^iV'2, (51) 

if in addition the free bosonic fields (p^ constitute potentials for axial vector currents 
composed out of the free Fermions ip^ 

S^0, = e.^J^ = ^„7^7'^„, A^O, a = 1,2. (52) 

The triple normal ordering in equation (^) is defined as ■.e'^'^: = e'^'^/ (^e'^^') for k 
being some constant. This is very advantageous in the calculation of commutation 
relations, since one can simply deal with ordinary operator relations instead of hav- 
ing to handle messy Wick contractions. We stress that in case the coupling constant 
A vanishes, that is when Cp reduces to £ff and the relations ( pOD correspond to 
two decoupled Dirac equations, the relation ( ^21) does not hold. 

In order to compute the factors of local commutativity 7^, as defined in (^, 
we need various (anti)-commutation relations. The fields ipai^) ^^^ complex free 
(Dirac) Fermions of masses rua and the fields 0„(x) are free Bosons, such that for 
a, (3 = 1,2 we trivially have 

[^ai^).My)] = [^aix),^f,{y)] = [0jx),$^(i/)] = {^Jx),V^^(i/)} = 0, (53) 

{^PM.^liy)} = S.p5ix'-y'). (54) 

The commutation relations involving mixed expressions of ip^ and 0^ are less obvi- 
ous and in fact it is crucial to note that these fields are not mutually local, that is 
[ipa{x),(pf3iy)] ¥" fo^^ space-like separations, i.e. (x — y)'^ < 0. Concretely we have 
the following equal time exchange relations for a, (3 = 1,2 

[^Ja;),0^(l/)] = V^5^pQ{x^-y^)^M. (55) 

^.(^)<^^(Z/) = $J(l/)^Jx)e2-(-i)'^^-^«(^^-^^), (56) 

-i^^{x)^!p{y) = vl/^(i/)^„(x)e-2-(-i)''^^i'^-^i.^«(-^-^^), (57) 

^a{x)^'^{y) = $^(|/)v^«(x)e2-^(-i)'^^"^®(^^-^'), (58) 

-^a{x)^(s{y) = ^/3(l/)^a(x)e~2"*^(-i)'^i'^-^i-i. (59) 
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We used here the Heavyside step function 6(x), defined as usual as 6(x > 0) = 1, 
0(x < 0) = and 6(0) = 1/2. One may convince oneself easily that ( ^5]) is compat- 
ible with (0) and that the remaining equations are straightforward consequences 
of (|53D-(^5D. Apart from this choice, which agrees with the one in ^^, one can 
also find in some places of the literature, e.g. |]31|, that in (|55| ) the 0-functions 
are replaced by e{x)/2 = Q{x) — <d{—x). This is of course also compatible with 
(^). However, an immediate consequence of our choice is that the Federbush fields 
"^ai^) are only mutually local if they are of the same type a, whereas when taking 
the ^-function instead, they are mutually local for all values of a and (3. The dif- 
ferent choices will of course lead to different factors of local commutativity 7 and 
will therefore alter the consistency equations (H)-®- Arguing on the properties of 
these equations we provide more reasoning for our choice below. 

A further important implication of the fact that ■?/'„ and 0^ are not mutually 
local is that the fields ^^q, are in different Borchers classes^] as the free Fermion. 
Thus, there is a chance for the existence of a nontrivial scattering matrix, which 



was indeed found in pB| , pT|] . In fact, we will now demonstrate that this S-matrix 
can be obtained as a limit of a more complex model, that is the homogeneous 
sine-Gordon (HSG) model. 

4.1 Federbush Models from HSG-models 

Ever since the equivalence between the massive Thirring- and the sine- Gordon 



model was demonstrated ||3^, there have been various identifications between dif- 
ferent types of models. In a similar spirit we also want to show now how a fermionic 
model is obtainable from a bosonic one, albeit in contrast to the above situation our 
fermionic scattering matrix will be constructed solely out of the asymptotic phases 
of a given scattering matrix 

hin 5.,(^) = e^l (60) 

Starting from a consistent solution to the crossing and unitarity relations 

Siti(^)S'di-0) = l and Sitie) = St{^7r - 6) , (61) 

one clearly has the constraint A^- = — A~^ = A^. This means that 

S,j = e^,+^3 (62) 

will also be a valid solution to the unitarity-crossing relations for S. Having no 
rapidity dependence there is no bound state bootstrap equation to be concerned 

^An equivalence class of complete, local field systems is referred to as a Borchers class. Its 
crucial property is that it characterizes completely the scattering matrix without having to resort 
to particular fields. (For more details see e.g. pOJ p. 104, however, this notion is of no further 
relevance for our concrete computations.) 
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about, such that ( |62D aheady constitutes a consistent scattering matrix. Concretely 



we will now show that when taking Sij{6) in (|60D to be the scattering matrix of the 
S'?7(3)3-HSG model, as found in [^, the resulting S-matrix Sij in ( |62D will be the 
one of the Federbush model at a particular value of the coupling constant. In fact 
it can be shown that this prescription leads to a much wider range of scattering 
matrices which can be directly associated to a Lagrangian of a Federbush model 
generalized in a Lie algebraic manner, see section 5. Let us recall now the scattering 
matrix of the S'f/(3)3-HSG model 



^5^(3)3 (^) 



/ (2). 

-(1)6 



V 



-(-2),e^ 
-l)ee-" 



-(1). 
(2). 



(2). 






(2). 



\ 



/ 



(63) 



We abbreviated {x)e = sinhi(6' + z7rs/3)/sinh |(6' — inx/3) and r = ±1/3. For 
the rows and columns we adopt here the ordering {1,1,2,2}. We also took the 
resonance parameters a of the HSG-model to be zero, since they will not play any 
role in our further considerations. Computing now the limit according to the above 
prescription we obtain 



lim[^f(^)^(^)^f(^)^ 



u 



S 



FB 



/ 


1 


1 


„— 27riA 


g2.*A 


1 


1 


^2niX 


^-2.^) 




27riA 


-27rJA 


1 


1 


\ 


-2mX 


27riA 


1 


1 



\ 



J 



■ (64) 



We found it convenient to relate the parameter r to the A in the Lagrangian density 
(^) as r = 1 — A. Then S^^ corresponds to the scattering matrix derived in 
p6| , |3T| , apart from the overall minus sign, which is due to the fact that we adopt the 
convention that the particles are ordered in opposite order in the in- and out-states, 
i.e. we include the statistics factor into the S-matrix. After having taken the limit 
(0), the crossing and unitarity equations also hold when we relax the constraint for 
T and allow it to take completely generic values different from 1/3. Thus, whenever 
the coupling constant A becomes an even integer the theory decouples into a system 
of two free complex Fermions. From a Lagrangian point of view we expect this kind 
of behaviour of course for vanishing A. 

Having specified the scattering matrix of the model, we are in the position to 
state directly from @ a representation for the FZ-algebra. The explicit version of 
(H) then reads 



Zl{e) = exp (-iX f d9':p^{e')-\ a\{e) 

z\{e) = exp (ix [ de'-.p^ie')-] a\{e) 



(65) 
(66) 
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Zl{e) = e^^Uxj^ de':p,{9')^a\{e) (67) 

Zl{e) = exp (-iX f de':p^{e')-] 4(^) (68) 

with 

PaiO) = aUe)a^{e) - aUe)a^ie) . (69) 

We will now specify more concretely various local operators of the Federbush model 
for which we want to compute the form factors explicitly by using the fermionic 
free field representation (|65D-(|68D. 

4.2 Form factors of some local operators 

We compute now explicitly the bosonic fields 4>a{x) by solving equation ( p^ and 
express them in terms of our general formula (|l^ 

/X — 

dx':^i^^: = xUx) with ^"(^, ^') = __I!-_^ . (70) 

A field closely related to 0Q,(a;), but whose origin is far less direct, is 

^'ai^)=XUx) With '^'(^'^') = -'i'M,-^') = ^|^^^^^. (71) 

It is this field which constitutes the analogue to the auxiliary field already used in 
the previous section. In view of the periodicity of the scattering matrix (^9]), we 
may restrict the range of A to A e (0, 2) /I. The special role of A = 1 was treated in 
more detail in ^^. Important for our purposes is the value A = 1/2 for which the 
operator ^^(x) reduces to uja{x) as defined in equation ([I^) . 



4.2.1 The order and disorder field 

In close relation to the free fermionic theory one may also introduce the analogue 
fields to the disorder and order fields in the Federbush model 

^'M = ■■ exp[fi^(x)] : and S^(a;) = :^„(x) $^(x): . (72) 



In [^, Lehmann and Stehr showed the remarkable fact that the operator $^(x), 
which is composed out of free Bosons, occurring in (^) can be viewed in two 
equivalent ways. On one hand it can be defined through a so-called triple ordered 
product and on the other hand by means of a conventional fermionic Wick ordered 
expression 

$^(x) = ; exp[-2v^zA0jx)]; =: exp[n^(a;)] : . (73) 
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Having again in mind to compute the factors of local commutativity 'jj^, as defined 



M • 



in (^, we need various equal time exchange relations. With the help of (p5|)- ([59|) 
we compute 

-^aix)^iy) = Sj(l/)^Ja:)e2-(-i)''^^"''«(^^-^^), (74) 

^'M^iiy) = ^',{y)K{x)e'-^^-'^'''-^^^^'-y'K (75) 

-^«(a;)S^(y) = sJ(y)^,(a;)e2"*^(-i)^(^"''®("'-^')-^i-^i-i®(^'-"')), (76) 

S^WSjd/) = Sj(y)S^(x)e^-(-i)'^^'^''. (77) 

Having obtained the relevant exchange relations we can read off the factors of local 
commutativity for the operators under consideration 

^*^ = -^f ^ = e2-(-i)''^^'^^ and ^l^ = -^f = e-'^'^-'^'^'-^ . (78) 

Note in particular, that for A -^ 1/2 we recover, as we expect, the values corre- 
sponding to the two complex free Fermions 

lim 7*" = 7^° = -1 and lim 7^- = 7^- = 1 . (79) 

A- — ^1/z A — >l/2 

Having assembled all the ingredients, let us now turn to the exphcit computation 
of the n-particle form factors related to the field $a(x). Since £ff respects the same 
symmetry as Cp, namely (|T3|), it is an immediate consequence that the only non- 
vanishing form factors of ^^(a;) have to involve an equal number of particles and 
anti-particles a and a. That means 

^^^,\aa...aaf3(,...(3f3^^^^ . . . ^fc, 9^+1 • • • , O^+l) = , « ^ ^, fc ^ /, «' = 1, 2, 1, 2 . (80) 

Turning now to the non-vanishing form factors, we compute by employing again 
Wick's theorem 



isin(7rA)e^^i2 ^ ^-^,22 



2 cosh ho 



FP '^^(01, ^2). (81) 
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Note, that in the contraction of an aa{0)- and a Z^(6')-operator there is no contri- 
bution from the exponential term inside the Z^{6), since it involves always particles 
of a different type than a, see (|65D-(|68D. Proceeding again in the same way as in 
the previous section, we obtain as closed expressions for the n-particle form factors 

^2n [Xl,X2 ■ ■ ■X2n-l,X2n) — [ — i-) J^ 2n [Xl, X2 ■ ■ ■ X2n-1, X2n) — 

*f^|nxll _ $^|nx22 _ 

F^n [Xi, X2 . . . X2n-1, X2n) = [-1) F^n [Xi, X2 ■ ■ ■ X2n-1, X2n) = 

^n^n-l sin"(vrA)a„(Xi . . . a;2„-i)^+^a„(x2 . . . X2n)^-^Bn,n ■ (82) 
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We may now convince ourselves, that the expressions for Fg^" " indeed satisfy the 
consistency equations (|1)-(|1). The first two equations are rather obvious to check 
and we will not report this computation here, but the verification of the kinematic 
residue equation @ deserves mentioning 

■D p*l[{"+l)xll/ - - - N _ 

-'^^§-^2n+2 l~3;, X, Xi, X2, . . . , X2n-1, a^2nj — 



''r+^ sin"+^(7rA)o-„+i(-x, xi . . . X2„_i)^+^o-„+i(a;, Xs . . . X2„) 2 ^ Res; 

x— ►x 

2n 

il>;^|nxll /_ 



k=l 



F^n (Xi, . . . X2n-1, a;2n) • (83) 



Recalling the definition of Bn,n of (^), we used 'ReSx^xBn+i,n+i = —x ^Bn,n and 

the value for 7^^ from (ffSl). Note the factor sin(7rA), which was originally found 
135[| , and which appears in our presentation in (^) relatively unmotivated, is 



m 



absolutely crucial for the validity of ([83|). 

Similarly we evaluate the matrix elements of E^ 



pE^|l{nxll).^ ^ ._. .„pS2-^|2(nx22).^ ^ 

-^2n+l 1^1, . . . ,y2n+lj — (, — ij -r2„+l (.f 1, • • • , f2n+l 

~E--^|l(nxIl)^^ ^ ^_^^^„~S^|2(nx22)^^ ^_ ^_(2i)" 



^;^">^^^\^„...,^2„.i) = hl)-FiST'''''i0u...,9.n.^) = ^sin"(.A) 



X 



Cr„(X2...X2n)^2 -TT /- - \ V^ J<l-J,¥k 



, . n (^2.-X2,)E TT^^'^^^^ • (84) 

a„(xi . . . X2„+i) 2 !<,<,•<„ fc (xk)^ n n(a;j + a^O 

- - j^k I 

However, the expressions of -P^n+i o^^Y satisfy the consistency equations (@)- 

(^ for A = 1/2. This reflects the fact that S^(x) is only a local operator for this 
value of A, see equation (|77D . Thus, the equations (H)-(|) "know" about the locality 
properties of the operator involved. 

As we already commented above, part of the operator content of the Federbush 
model reduces to the one of the complex fermionic theory. We may check explicitly 
that the same limit is respected by the form factors 

lim Fi^l"^"-(ei,...,M = F,':r'"""(^i,---,M, (85) 

lim ^2^.^1^""^"^(^l,---,^2„+l) = F^4f"'^"^(^i,...,^2n+i). (86) 

A — *l/z 

Note, however, that since the real fermionic theory can not be obtained directly 
from the Federbush model, see also section 3.1.1., we also do not recover, as we 
expect, the same expressions for the form factors when the particles are taken to 
be self-conjugate. 
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4.2.2 The Federbush fields 

Let us now compute the form factors of some fields which occur explicitly in the 
Federbush model. From the expressions of the previous section the form factors for 
the Federbush fields follow easily 

(87) 



2n+l If^l) • 


•,^2n+l) = 




■ ,d2n) Ui{92n+l) , 


2n+l ic/i, . 


■,02n+l) = 


= V7r^2n i^U- 


• ,d2n) U2{62n+l) ■ 



Recall the definition of the Weyl spinors Ua{0) from equation (fi^ . It is clear that 
for each component these fields satisfy the form factor consistency equations. As 
already mentioned in section 3.1.2. and as is quite common in the literature, e.g. 
|]37| , ^ , third reference in p etc., one may verify that various equations which hold 
for the operators are also satisfied by the related form factors. However, one should 
be aware that such relations also hold for the matrix elements F, which do not yet 
satisfy the consistency equations (0)-(Q). Hence, the only conclusion one may draw 
from such comparisons is a relative consistency amongst the solutions obtained. 
Such arguments do not serve as a stringent identification of the operators, albeit 
they give an indication. We illustrate this statement with the following simple 
computation. Let us take the fields as defined in (|TDp and evaluate directly by 
Wick contracting 

Ff'^'^l""(^,^) = v^F/^l""(^,^) = -i7rim„cosh^, « = 1,2, (89) 
p9i^J^»^0^0^ = -y^F/°l""(^,^) = z7rim„sinh^^, a = 1,2. (90) 



This confirms precisely the conservation equations (|52D on the level of the matrix 
elements. However, it is also easy to see that the expressions (^) and (^) are 
not yet solutions of the form factor consistency equations (0)-(|^). In principle, 
these equations together with the Dirac equation already ensure that the "^a are 
solutions of the equations of motion (|5(]|). Nonetheless, it is instructive to verify 
(pop explicitly. Using still the representation (|70D, we compute 

Ft'^'"'''^''\eue2,e,) = z7clm2u,{d^ + 62-63) (91) 

^,7^9^*11221.. n n\ ^^^^ Km2Ui{6i3 + 6i) + m2Ui{6i3 + 6i) + miUi{63)) 

^3 [61,62,63) = -^7 (92) 

cosh 26^12 

f:^^''\6„62,63) = ^^u,i63) . (93) 

cosh 7^612 



Assembling these expressions, we confirm directly the validity of (^OD at the level of 



the three particle matrix elements. We expect of course this property also to holds 
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for higher orders. It is easy to check that (]91|)-(]93|) do not constitute solutions of the 
equations (H)-®, in particular (|87|) does not reduce to (|93|) . Thus on one hand we 
see that formal operator equations do not serve as a conclusive means of operator 
identification and we therefore need alternative arguments such as the ultraviolet 
limit in section 5 etc. On the other hand this underlines further the need for the 
introduction of the field f2^fx). 



4.2.3 The energy momentum tensor 



The energy-momentum tensor for the Federbush model has been computed in ||39|| . 
Its evaluation involved a small subtlety, since the one obtained directly from the La- 
grangian does not lead to the correct Poincare generators, such as (^Tj). This could 
be fixed in the usual way by exploiting the ambiguity in the definition. Essential 
for our purposes is once again the trace, which is 



T'^ = 2imi:^i^i: +2im2:^2^2: 



(94) 



Using the representation (|5TD for the Federbush fields, we compute the only non- 
vanishing form factor for T^ to 



"T-t Li 



\9,~9) 



-.T^Jaa, 



-2'Kim'i, sinh 



9-9 



(95) 



This means the function is the same as the one for the complex free Fermion. 



4.3 Momentum space cluster properties 

As a consequence of Weinberg's power counting theorem one has also a further 
property of form factors which involves the structure of the operators themselves, 
namely the momentum space cluster property, see e.g. [|l| some reasoning on this. It 
serves on one hand as a consistency check for possible solutions of (0)-(§) and on the 
other as a construction principle for new solutions, e.g. |36|. It states that whenever 



some of the rapidities, say k, are shifted to plus or minus infinity, the n-particle 
form factor related to a local operator O factorizes into a k and an {n — K)-particle 
form factor which are possibly related to different types of operators O' and O" . 
Introducing the translation operator T^ which acts on a function of n variables as 



Ta fi^li- ..,9a,...,9r, 



f{0 



li 



,9a+^,...,9r, 



(96) 



and the operators 






i™n^".. 



p=a 



and 






lim Y\T 



±1? 

2p ' 



p=a 



(97) 
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the statement of inomentuiii space cluster decomposition reads 



Ta.KTa.KFn (6*1 •• • 6^^ 



"2(K-a+l)( 



72a- 



-1 • • • d2K)Fn 



O'' 



-2(K-a+l)(^l • ■ ■ ^2a-2, ^2k+1 ■ ■ ■ On) ■ 

(98) 

Of course, we could have defined the product of T^^T^^ to be just one operator, 
but it will be convenient for us to distinguish the shifts in even and odd positions 
of the particles. Let us now see the effect of the action of these operators on the 
various functions which build up our form factor solutions, see (^7|)' (H) ^^"^ (HI)- 
We compute 



■±rr± 






and 



1^+5 



(Jn\X\ . . . X2n-1 
amix2 . . . X2m)^''^ 



1,K 1,C "'™ '^ ^K,<;i^n-K,; 



~ e 



±M{K-C)±-d^^^ 



an(xi 



(Jr 



X2„_l)^+^ 



{x2...X2mV 2 



&«-") 



■k-<: 






-^Ol^+^Oi + O 



2? J 

t(^2K + l---^2n-l) 1 
-?(^25+2---^2m) J 



(99) 



(100) 



In order not to overload our symbols, we have slightly abused here the notation. 
Whereas in (^) the Xj, x^-dependence of Bn,m always start at i,j = 1, in ( |100|) 
the dependence of Bn-K,m~q for the minus shift is the same as in the corresponding 
factor for the symmetric polynomials. Besides the explicit functional dependence 
on the r.h.s. of (^) and ( |100| ) it is instructive to consider at first the leading order 
behaviour 






anjXl . . .X2n-l) 
(^m[X2 ■ ■ ■ X2m) 



A+i 



J(C-'^)\ 
Bn,m ~ ^ i5{C-f=)| 



-n+JTi — A 



(101) 



From this we see directly that in general the final expression will tend to zero, unless 
( = K, \( — K,\ = 2X or \( — K ± 2\ = 2A, by noting that our solutions only allow 
\n — m\ = 0, 1. So, let us now collect the functional dependences in equations (|99|) 
and ( |100| ) and see how our form factor solutions combine under clustering to new 
form factors. We compute 






<J>Q|nxaa 
2n 



p<S>^\KXaa,Q 






2n 



-'l,K+l^l,K-^2n 



-'l,K-'l,re+l-' 2n 
'J-+ '1'+ rp'^a\(nxaa)a ■p^i^\KXaa / 



2k 

2k \^l 

<Ja \{KXaa.)a / / 
2k+1 V- 

-^<j a\{KXaa)ci / / 
"2k+1 V- 

2k 



.e,^F: 



Fr 



72KK2(n-K) y 

^2KJ^2(n-K) I 



^2k+1 • • • ^2nj 
^2k+1 • • • ^2n) 



'1 



(102) 
(103) 



1 • • • ^2K+l)F-'l^:)~_7 ^ (^2K+2 . . . ^2n), (104) 

1...6'2K+lj-^2(n-K)-l [ty2n+2---t>2n), (105) 

• • ^2K)F2TJ-KT+f ''"'"(^2'^+l • • • ^2n+l) , (106) 
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q— n — T^(Ta\{nxaa)a 


rsj 


pii.^\KY.aa,Q 


/I X pO-cl[(n-K)xaa]a/^ 
^2K)J^2(n-K)+l l^2«+l • • • 


d2n+l) , 


(107) 


/J-+ q-+ po-c 1 (n X aa)a 
^l,K^l,K+l-^2n+l 


r*^ 


^2k+1 1^1 


n \ T7Mal{"-K)xaa/^ 
• • •y2K+lj-^2(n-K) Id'2k+2 


■ ■ ■ ^2n+l), 


(108) 


q-- q — pCTa (nxaa)a 
^l,K^l,K+l-^2n+l 


rsj 


^2k+1 1^1 


• • •y2K+lJ-r2(„_K) lf72K+2 


• • • (^2n+l), 


(109) 


q-+ q-+ po"a (nxaa)a 
^l,K^l,K-^2n+l 


r*^ 


^2k 1^1 • • 


/] \ T70-ii|[(n-K)xaa]a//i 
f^2KJ^2{n-«)+l l^2«+l • • • 


02n+l) , 


(110) 


^l,K^l,K-^2n+l 


r*^ 


^2k 1^1 • • 


n \-rp(Ta\[{n-K)xaa\a,i;. 
f^2K)^2(n~K)+l {^2k+1 ■ ■ ■ 


^2n+l) , 


(111) 


'7-+ 'r+ pO"a|{"Xoa)a 
^l,K+l^l,K-'^2n+l 


rsj 


i-l(Ta {KXOa)a //J 


• • • ^2K+lj^2{n-K) lt^2K+2 


• • • ^2n+l), 


(112) 


^l,K+l^l,K-^2n+l 


rsj 


j-i(Ta|(KXaa)a/^ 

^2k+1 1^1 


n N T7/ial{"-K)xaa!//i 
. . . 6'2K+lj-^2{n-K) \.^2k+2 


• • • ^2n+l)- 


(113) 



Thus, omitting the shift operators we have formally the following decomposition of 
the operators 



$' 



a a 



CTr 






f^c 



fJ'ct X fJ'a 



(114) 



together with the equations for a ^ a. This means the stated operator content 
closes consistently under the action of the cluster decomposition operators. 

5 Lie algebraically coupled Federbush models 



The Federbush model as investigated in the previous section only contains two types 
of particles. In this section we propose a new Lagrangian, which admits a much 
larger particle content. The theories are not yet as complex as the HSG-models, 
but they can also be obtained from them in a certain limit such that they will 
always constitute a benchmark for these class theories. Form factors related to 
these models may be computed similarly as in the previous section. 

Let us consider i x £-real (Majorana) free Fermions ■?/'^ (x), now labeled by two 
quantum numbers l<a<i, l<j<i and described by the Dirac Lagrangian 
density £ff- We perturb this system with a bilinear term in the vector currents 
J^ . = ^ ■^'^\l) ■ 



e i it 



(115) 



a=l j=l 



--lj,k=l 



and denote the new fields in Cqf by ^aj- Furthermore, we introduced £^ x i"^ 
dimensional coupling constant dependent matrix A^^, whose further properties we 
leave unspecified at this stage. As in the usual Federbush model, the effect of the 
presence of the Levi-Civita pseudotensor e is that the theory described by £cf is 
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not parity invariant. Thus £cf may be viewed as a system of coupled Federbush 
models p5| . 

The formal equations of motion associated to Cqf are easily derived as 

e i 
{lYd, - m,,)vl/,, = 7re,.7^ E E ^a'^/.^fc^a, . (116) 

b=l fc=l 

The solutions to these equations can be constructed in close analogy to the ones of 
the Federbush model. The fields 



^a,j = ■ exp(v/^Z Yl Yl Kl<f>b,ky-^a,j = ^iji^aj (117) 



solve the equations of motion ( |116| ) with the additional assumption that the bosonic 
fields 0^j- constitute potentials for axial vector currents 

^9,0,,,. = e,,J:^^ = ^,,,7^7'V'a,,, Kl ^0, yb,k. (118) 

As in the previous section, we used here once again the triple normal ordering in 
equation ( |117| ). It needs further computations, similar to the ones for the Federbush 
model, to make it rigorous that also in this context the triple ordering can be 
associated to a standard Wick normal ordering. Nonetheless, it appears natural 
to expect that this can be generalized analogously and we take this here as an 
assumption. 

Accepting this, we can now compute various equal time exchange relations with 
l<a,b<i,l<j,k<i 

[(f>aA^)Ab,k{y)] = [^aA^),^bAy)] = o (119) 

bPa,jix)Ab,kiy)] = \^Sa,bSj,kQ{x^ -y^)'iPa,jix) (121) 

i^aA^Hkiy) = Kk{y)4^aA^)e''''''^^'"^^'-''^ (122) 

-ij,A^)^,Ay) = ^6,.(2/)^.,,(a:)e-^-^-®(-^-^') (123) 

-^aA^)'^b,k{y) = *6,fe(2/)^a,i(x)e-^-^-. (125) 

The equations ( |119[ ) and ( |120| ) are again clear since ^^ • and 0^ • are free Fermions 
and Bosons, respectively. Equation (|121|) is compatible with ( |118| ) and the re- 



maining equations are simply consequences of ( |119| )-( p!2lD . With the help of these 
equations we compute directly the scattering matrix. We will be slightly casual 

22 



here about complete rigour and do not worry with test functions and smeared out 
operators. Noting that 0„j |0) = 0, we obtain from 

hm ^a,,^,,, |0) = Sitij^^^ij,^, |0) , hm ^,,^,,, |0) = ^^^..fcV'a,, |0) (126) 

the S-matrix 

Sit = iSitr'Sit = -e-^- . (127) 

Let us now see whether ( p.27| ) is consistent in the usual sense, i.e. that it passes all 
the tests of consistency or if the latter put some constraints on the possible values 
for the coupling constant dependent matrix A;^^,. We demand the usual crossing and 
unitarity relations (|6lD , which means we should have 

Ait = -A^i + 2Z and A^^ = A| + 2Z . (128) 

We will now provide some concrete solutions to ( |128| ) and therefore ( |127| ). 



5.1 HSG-type solutions 

Let us take 

Ait = 2Kte,J,,Kj (129) 

where K denotes the Cartan matrix of SU{N) and / the incidence matrix of a 
simply laced Lie algebra, which we refer to as g. The Xab sue i"^ coupling constants, 
which are, however, not entirely independent of each other. For instance we assume 
Aafe = Afea. Furthermore, we characterise the anti-particle exclusively by the first 
quantum number, i.e. (a, i) = (a, i), where the particle a may be constructed from 
a by the automorphism which leaves the associated Dynkin diagram invariant. In 
the case of SU{N), we simply have a = i+1 — a. It is clear that ( p.29| ) satisfies the 
first relation in ( |128| ), whereas the second relation introduces further constraints 



on the A's. To be more concrete we specify now ( |1271 ) for some special choices of 
A^ and the Lie algebra g. 

5.2 The Federbush model 

Considering now the case SU{3)3 with An = A22 = — 2Ai2 = — 2A21 = A, we obtain 
the scattering matrix of the Federbush model S^^ as defined in (|64D , where we 
now used the ordering {(1, 1), (2, 1), (1, 2), (2, 2)}. In comparison with the previous 
section, one should notice, that we have now realised this model in terms real 
Fermions rather than complex ones. 
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5.3 % 



To illustrate the formulae (|128|) and ( |127| ) a bit more, let us consider a sliglily 



more complex model, namely g^. When specifying the quantities in ( |129| ) to these 
algebras, we obtain 



S^3 



/ 


^2-Ki\eijIij 


^2mX'EijIij 


1 


^-2-Ki\' Eijhj 


^-2niXsijIij 


\ 




^-2TviX'eijIij 


g27rJA Bijhj 


1 


g-27riA e.jhj 


^2-ITi\'£,jUj 






1 


1 


1 


1 


1 






^2m\'eijlij 


g-27riA Eijlij 


1 


g27riA Eijlij 


^-2-Ki\'£ijIij 




\ 


^-2lTi\£ijIij 


^-2m\'eijUj 


1 


^2mX'EijIij 


Q^TTiXSijIij 


) 



(130) 



where the rows and columns are ordered as 1,2,... 5. In this case, we have three 
independent coupling constants A, A' and A". 

5.4 The HSG-limit 

From the specific example in (|6^ , we expect that the HSG-models are in general 
closely related to ( |127|) with ( |129| ). Indeed, taking \ab = 1 for 1 < a, 6 < £, we 
obtain 

Sll = -e^-i'^^hkK:^ ^ (131) 



which clearly satisfies the first relation in (|128|) , whereas for the second relation, 
we simply have to recall the well-known fact that {Kg^,j^C)ah = min(a, 6) — ab/N. 
Comparing with the expression 

hm K'"'^(^)^2"'^(-e)] = e2-^(^"'"')-^^ . (132) 

we note that these solutions coincide. This means when one eventually solves the 
HSG-models, one can always take the limit to the corresponding quantities of £cf 
for consistency checks. 



6 The ultraviolet limit 

When having found a solution to the form factor consistency equations, with the 
factor of local commutativity and the scattering matrix as the only input, one 
normally does not know which operator this particular solution corresponds to. Of 
course in the present situation we are in a better position, since we are already 
working with an explicit representation for the operators. Nonetheless, in section 
4.2.2. we saw that even this can still lead to wrong assignments and it is desirable 
to have more information. By calling the operators $a,/i, and S^,o", disorder 
and order operators, respectively, we have already borrowed the terminology from 
the underlying conformal field theory. In order to make this correspondence more 
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manifest one may carry out explicitly the ultraviolet limit. The ultraviolet Virasoro 
central charge of the theory itself can be computed from the knowledge of the form 



factors of the trace of the energy-momentum tensor |23] by means of the expansion 



^.k"'P""i ■■■! (EL.™. cosh..)* " **•■■■•" 



— oo 



(133) 
In a similar way one may compute the scaling dimension of the operator O from 



the knowledge of its n-particle form factors |24 



- /n\ Z^ Z^ 



dui . . . duji 

J ■■■ I z;;^;:^, 

n=l a-, ...u„ 

f 1 fn_oo 



2 {O) /-^^^ 7 • • • J ^,(2^)n [Y^Um,^ cosh^O^ 



X 



irniM,.../.„(^^^ __^^j (Ff^i-'^"(^i,...,^„))* . (134) 



In general the expressions ( |133|) and ( |134|) yield the difference between the cor- 



responding infrared and ultraviolet values, but we assumed here already that the 
theory is purely massive such that the infrared contribution vanishes. Let us now 
evaluate these formulae. 

6.1 The complex Fermion 

Since for the free Fermion one only has to sum up to the two particle contribution, 
the infinite sum (|133|) and ([L34|) terminate and can be evaluated even analytically. 



For the case A^ = 2 we obtain 

1 



,. = 2 and 1^^, = ^-, = —. (135) 

lo 



The scaling dimensions of a a and a a, which are expected to coincide with (|135|) , 
can not be computed from ( [134| ), since it involves an odd number of particles. 

6.2 The Federbush model 

We may proceed similarly for the Federbush model. In the ultraviolet limit it 
obviously corresponds to two complex free Fermions and we there expect to obtain 

Cu, = 2 . (136) 

Indeed using ([95|), the computation is identical to the one carried out in the previous 
section. Note, that this value of 2 coincides with the ultraviolet central charge of 
the 5't/(3)3-HSG model. This is however also not entirely surprising by recalling 
the identification (0). The corresponding thermodynamic Bethe ansatz equations 
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will be identical to for free Fermions. More striking is the result of the evaluation 
of (|134D, which yields with (H) and (|8l|) 



A*l = A*l = ^ . (137) 

1/2 1/2 

Note, that Z^^ = A**^ = 1/16, which is once again the limit to the complex 
free Fermion. Yet more support for the relation between the S'f/(3)3-HSG model 
and the Federbush model comes from the analysis of A = 2/3, which corresponds 
to the S'[/(3)3-HSG value r = 1/3 (see (|63|) ). In that case we obtain from (|137|) 



».2/3 ,2/3 



A*"^ = /x[^ = 1/9. We now compare with the general formula for the scaling 
dimensions of the S't/(3)3-HSG model 

A(A..).(^i4±M_(^, (138) 

1/ o 

where A is a highest weight vector of level smaller or equal 3, w the corresponding 
lower weights and p the Weyl vector. We are specially interested in the field cor- 
responding to A(Ai, Ai) with Ai being a fundamental weight, since this field was 
previously observed [Q to correspond to the disorder operator. Indeed, we find 
that 

»2/3 »2/3 

A(A„AO=A*i =A*i . (139) 

Thus precisely at the value of the coupling constant of the Federbush model at 
which the S'f/(3)3-HSG S-matrix reduces in the limit (Q) to the S^^, the operator 
content of the two models overlaps. 

7 Conclusions 

We computed explicitly the form factors of the complex free Fermion and the Feder- 
bush model related to various operators. On one hand we carried out this task by 
representing explicitly the field content as well as the particle creation operator in 
terms of fermionic Fock operators and computed thereafter directly the correspond- 
ing matrix elements. On the other hand we verified that these expressions satisfy 
the form factor consistency equations only when the operators under consideration 
are mutually local. This can already be seen for the free Fermion, for which we 
could have also computed the matrix element of the field $^(a;). In that context 
one observes that only for A = 1/2 this function solves the consistency equations 
(0)-(^. We observed a similar phenomenon in the Federbush model. Whereas the 
matrix elements of the field ^^(a;) can be computed in a closed form for generic 
values of A, they become only meaningful form factors for A = 1/2, that is when 
the field becomes local. It turned out to be crucial that the consistency equations 
contain the factor of local commutativity as defined in (^). It is important to note 
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that this factor is related to the equal time exchange relation between the operator 
and the field associated to the particle next to it in the multi-particle state. 

Our solutions turned out to decompose consistently under the momentum space 
cluster property. 

Further support for the identification of the solutions was given by an analysis 
of the ultraviolet limit. 

We demonstrated how the scattering matrix of the Federbush model can be 
obtained as a limit of the S'f/(3)3-HSG scattering matrix. This "correspondence" 
also holds for the central charge, which equals 2 in both cases, and the scaling 
dimension of the disorder operator at a certain value of the coupling constant. We 
proposed a Lie algebraic generalization of the Federbush models, which on the other 
hand can be obtained in a certain limit of the homogeneous sine-Gordon models. 

We expect that the construction of form factors by means of free fermionic Fock 
fields can be extended to other models. 
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